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Area preserving maps provide the simplest and most accurate means to visualize and quan- 
tify the behavior of nonlinear systems. Convenience of the mapping equations of motion for 
investigation of transition to chaotic behavior in dynamics of classical atom in microwave 
field, transition to nonchaotic behavior in randomly driven systems and induced quantum 
dynamics of simple and multilevel systems is demonstrated. 



1 Introduction 

A common method of displaying the dynamics is through a Poincare section. The Poincare 
section is a device invented by Henri Poincare as a means of simplifying phase space diagrams 
of complicated systems. It is constructed by viewing the phase space diagram stroboscopi- 
cally in such a way that the motion is observed periodically. A dynamical system whose phase 
space is tree-dimensional may be converted through the Poincare section to a two-dimensional 
mapping. Sometimes such mappings may be further simplified to one-dimensional. Because 
of their relative simplicity one- and two-dimensional maps provide several advantages over 
the differential equations. They allow for simple reveal of many characteristics of chaotic 
behavior, such as sensitivity to initial conditions, illustration the mechanisms of bifurcation 
and so on. 

Area preserving maps provide the simplest and most accurate means to visualize and 
quantify the behavior of conservative systems. Such maps may be iterated on even the simple 
computer. For analysis of the non-conservative systems sometimes the non-area-preserving 
maps may be introduced. 

Here we demonstrate examples of investigation of transition to chaotic behavior in the 
dynamics of classical atom in microwave field by area-preserving maps, transition to non- 
chaotic behavior in randomly driven systems by the non-area-preserving maps and induced 
quantum dynamics of simple and multilevel systems represented by the appropriate quantum 
maps. 
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2 Kepler maps 



Consider dynamics of the classical hydrogenic atom in the monochromatic field. The direct 
way of coupling the electromagnetic field to the electron Hamiltonian is through the A • P 
interaction, where A is the vector potential of the field and P is the generalized momentum of 
the electron. The Hamiltonian of the hydrogen atom in a linearly polarized field F cos(u;i+i?), 
with F, uj and •& being the field strength amplitude, field frequency and phase, respectively, 
in atomic units is 

H = - (p — — sm(cot + J --. (2.1) 

Electron energy change due to interaction with the external field follows from the Hamilto- 
nian equations of motion [1] 

E = -i-¥coa(ut + &). (2.2) 

Using parametric equations of motion in the Coulomb field we may calculate the change of 
the electron's energy in the classical perturbation theory approximation. 

Measuring the time of the field action in the field periods one may introduce the scale 
transformation where the scaled field strength and the scaled energy are F s = F/oj 4 ^ and 
E s = E/u 2 / 3 , respectively. However, it is convenient [2, 3] to introduce the positive scaled 
energy e = —2E S and the relative field strength F = Fuq = F s /sq, with n being the initial 
effective principle quantum number, n = (— 2E o y 1 ^ 2 . The threshold values of the relative 
field strength F for the ionization onset depends weaker upon the initial effective principle 
quantum number n and the relative frequency of the field s = urn^ than the scaled field 
strength F s . 

We restrict our subsequent consideration to the one-dimensional model, which corre- 
sponds to the states of low orbital quantum numbers / <C n and is widely used in theoretical 
analysis [4-8]. Integration of Eq. (2.2) for motion between two subsequent passages at the 
aphelion (where x — and there is no energy exchange between the field and the atom) 
results to the map (see [9-11] for details) 



Ej + i = Sj — 7iF elh (sj+i) sint?j, 

= fij + 2irejl{ 2 - nF elr] (e j+1 ) cos i? 



3/2 " _ ip /_ \ o ( 2 - 3 ) 



where 



h(e j+1 ) = -^3'(s j+1 ). (2.4) 



£ ii s ^ +1 



Here s = e 3 / 2 = 2E) 3 ^ 2 = uj/VL is the relative frequency of the field, i.e., the ration 
of the field frequency u to the Kepler orbital frequency Vl = (—2E) 3 ^ 2 , and J' s (z) is the 
derivative of the Anger function with respect to the argument z. The function r\ (£j+i) may 
be obtained from requirement of area-preserving of the map (2.3) 

This requirement yields 

^.f^ii). (26) 

ae j+1 
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The derivative of the Anger function 

1 u 

J' s (s) = — / sin [s (x — sin re)] sinxdx (2.7) 

7T J 


is a very simple analytical function which may be approximated quite well by some combi- 
nation [10] of expansion in powers of s 

J >W = 2 lt-*) aiR ™' - 3 - 1 (2 ' 8) 

and of the asymptotic form 

,, . . b a sin7rs _ _ _ 

S 2 / 3 5S 4 / 3 47TS 2 

where 

2i/3 2 2 / 3 
a = — — — ~ 0.4473, b = — -—- ~ 0.41085. (2.10) 
3 2 / 3 r(2/3) ' 3^(1/3) V ; 

The map (2.3) is the general mapping form of the classical equations of motion for the one- 
dimensional hydrogen atom in a microwave field derived in the classical perturbation theory 
approximation. Some analytical and numerical analysis of this map has been done in Refs. 
[9-11]. Here we analyze different special cases of the map (2.3). 



2.1 High frequency limit 

For the relatively high frequencies of the field, s » 1 (s > 2), theoretical analysis of the 
classical dynamics of the one-dimensional hydrogen atom in a microwave field is relatively 
simple. That is why, the energy changes of the electron, (Ej + i — Ej) and (£j+i — £j), do not 
depend on the initial energy Ej and relative frequency s ^> 1. Indeed, using the asymptotic 
form of the derivative of the Anger function, J' s (s) = b/s 2 ^ 3 , we have h (£j+i) = Ab = const., 
r] (£j+i) = and, consequently, the following map 




0j + 2tt£- + 3 { 



3/2 (2.H) 



Note, that scaled classical dynamics according to maps (2.3) and (2.11) depends only on 
single combination of the field parameters, i.e., on the scaled field strength F s = F el = 
F/cu 4 / 3 . 

By the standard [12, 13] linearization procedure, Ej = eo + Aej, in the vicinity of the 
integer relative frequency (resonance), s — £ o 3 ^ 2 — m with m integer, the map (2.11) may 
be transformed to the standard (Chirikov) map 

Here Ij = -3irAej/e 5 /2 and K = l2n 2 bF /^. 
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From condition of the onset of classical chaos for the standard map, K > K c ~ 0.9816 
[12-15], we may, therefore, estimate the threshold field strength for chaotization of dynamics 
and ionization of the atom in the high frequency field 

F C = KJ (12tt 2 6so /3 ) ~ 0.02s 1/3 . (2.13) 

Some times [6] one writes the map (2.11) for a variable N = — l/2n 2 uj, the change of which 
gives the number of absorbed photons, 

N j+1 = N J +2n(F/u J ^)smrJ v 
= + 2nu (-2iuN j+i y 3/2 . 

We see that for such variables dynamics of the system depends on two parameters: on the 
quantum scaled field strength F q = F/uj 5 ^ 3 [2, 3] and on the field frequency u. Map (2.14) 
is, therefore, not the most convenient one for analysis of the classical dynamics. 

In general there are, however, no essential difficulties in theoretical analysis of classical 
nonlinear dynamics of the highly excited hydrogen atom in the microwave field of relative 
frequency s = oon^ > 0.5 when the field strength is lower or comparable with the threshold 
field strength for the onset of classical chaos, i.e., if the microwave field is considerably weaker 
than the characteristic Coulomb field. In such a case, energy change of the electron during 
the period of intrinsic motion is relatively small and application of the classical perturbation 
theory for derivation of the Kepler map (2.3) is sufficiently correct. Further analysis of 
transition to chaotic behavior and of the ionization process may be based on the map (2.3) 
and for s — 0.3 -j- 1.5 results in the impressive agreement between measured ionization 
curves and those obtained from the map (2.3). Even analytical estimation of the threshold 
field strengths based on this map is rather proper [9-11]. 

Sufficiently more complicated is analysis of transition to stochastic motion and of ioniza- 
tion process in the region of low relative frequencies, s < 0.3. 



2.2 Low frequency limit 

For the low relative frequencies of the microwave field, s < 1, the map (2.3) may be simplified 
as well. Using expansion of the function J' s (s) in powers of s, J' s (s) ~ s/2, for s < 1 we 
have according to Eqs. (2.4) and (2.6) 



h(e j+1 ) = 2/ef+ 1 
V = -5/4+1- 



(2.15) 



Consequently map (2.3) transforms to the form 



= A, + 27T/4 2 , + 57TF, (Eg/^,) cos d r 



This map is a little bit more complicated than map (2.11) for high frequencies, however, it 
may easily be analyzed as numerically as well as analytically. Note first of all, that energy 
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change of the electron during the period of intrinsic motion (after one step of iteration), 
\£j+i — £j \ , is considerably smaller than the binding energy of the electron Sj ~ e if the 
field strength is lower or comparable with the threshold field strength, i.e., 2nF ^£o/ £ j+i) — 

2vrF e 1/2 < e , or 27rF s < 1 if F < - °- 13 and s o < 1- This indicates that the 
map (2.16) is probably suitable for description of dynamics even in the low frequency region 
where the field is relatively strong. 

2.2.1 Adiabatic ionization 

For low frequencies, 2ixs = 2n/e 3 ^ 2 <C 1, according to the second equation of map (2.16) the 
change of the angle d after one step of iteration is small. As it was noticed above, the energy 
change is relatively small too. Therefore, we may transform difference equations (2.16) to 
differential equations of the form 

= -H^sintf, 

27r £ . 57r £ 2F q ( 2 -17) 

= + ^77^COS$. 

Dividing second equation of the system (2.17) by the first one we obtain one differential 
equation 

d(cos$) e 5cos-# 

Analytical solution of Eq. (2.18) with the initial condition e = Eq when $ = $0 is 

costf = z 5 cos<d -2z A (l- z) /F , z =^e/e Q . (2.19) 

Eq. (2.19) describes motion of the system in e and d variables, i.e., represents functional 
dependence between two dynamical variables. For relatively low values of F , i.e., for F < 
^z 4 = I , the right-hand side of Eq. (2.18) is positive for all phases ■&. Therefore, cos-# 
and e decrease and increase simultaneously and, according to Eq. (2.19), there is a motion 
in all interval [0, 27r] of the angle 1?. For F Q > |z 4 , however, the increase of the angle 1? in 
the interval 4- n turns at t? ~ n into the decrease. This results in fast decrease of e and 
to ionization process. It is easy to understand from analysis of Eq. (2.19) that the minimal 
value of F for such a motion (resulting in ionization) corresponds to $0 = and 1? = ir. This 
value of F is very close to the maximal value of F Q resulting to the motion in all interval 
[0, 2n] of i.e., the maximum of the expression 

F = 2z(l-z)/(l + z 5 ) . (2.20) 

This maximum is at z = zo, where zo is a solution of the equation z 5 + 5z — 4 = 0, being 
z ~ 0.75193. The critical value of the relative field strength, therefore, is F ° = 2z^/5 = 
0.1279 which is only 1 lower the adiabatic ionization threshold F si = 2 10 / (3tt) 4 = 0.1298. 

2.2.2 Chaotic ionization 

For higher relative frequencies, so > 0.1, ionization process is due to chaotic dynamics of 
the highly excited electron of the hydrogenic atom in a microwave field. There are different 
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criterions for estimation of the parameters when dynamics of the nonlinear system becomes 
chaotic. For analysis of transition to chaotic behavior of the motion described by maps (2.3), 
(2.11) and (2.16) the most proper, to is the criterion related with the chaotization of the 
phases [13] 

K = max — £± - 1 > 1. (2.21) 



to; 

Here max means the maximum with respect to the phase i?j and variation of the phase flj+i 
with respect to the phase dj means the full variation including dependence of on -&j 
through the variable e j+1 in Eqs. (2.3), (2.11) and (2.16). 

Applying criterion (2.21) to the general map (2.3)-(2.4) we obtain the threshold field 
strength 

£ 7/2 

If e ~ e Eq. (2.22) yields the result 

F c =(l27r 2 S J' s ( S )) _1 (2.23) 

which for s ^> 1 coincides with Eq. (2.13). 

For more precise evaluation of the critical field strengths we should take into account the 
change (increase) of the electron's energy due to the influence of the electromagnetic field 
[11]. For higher relative frequency s or lower scaled energy Sj the threshold ionization field is 
lower. Therefore, if the scaled energy Sj decreases in a result of relatively regular dynamics 
in not very strong microwave field, then it is sufficient the lower field strength for transition 
to the chaotic dynamics. 



3 Transition to nonchaotic dynamics and synchroniza- 
tion in randomly driven systems 

When an ensemble of bounded in a fixed external potential particles with different initial 
conditions are driven by an identical sequence of random forces, the ensemble of trajectories 
may become identical at long times, i.e. synchronization of the identical systems by common 
noise may be observed. Fahy and Hamann [16] considered a particle of mass m moving 
according to Newton's equations in a potential V(x), except that at regular time intervals 
r the particle is stopped and its velocity is reset to random value chosen from a Maxwell 
distribution with temperature T . It should be stressed that for every particle of the ensemble 
it was given an identical, randomly chosen velocity at the start of each step of time length r. 
This motion is in many respects similar to Brownian motion of the particles at a temperature 
T. However, if the time interval r between stops is lower than a threshold value r c , the final 
trajectories of the particles are independent on the initial conditions; all trajectories become 
point by point identical in time. Although the trajectory is highly erratic and random, the 
system is not chaotic. 

The similar effect may also be observed in a more general and realistic (from the physical 
point of view) case, i.e., when mixing at time intervals r some part a of the old velocity 
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\ old with random velocity v mra to get a new starting velocity v r 



av 



old 



Here a 



threshold value r c depends on a. 

Let us consider a particle of mass m moving in a one-dimensional potential V(x) which 
confines particles to a finite region. At a time intervals Tj the particle is partially stopped 
and its velocity is reset to a new starting velocity Vi 
particle moves according to Newton's equations 



avf d + vl an . Between the stops the 



dx 
~di 



d 2 x 1 dV 
dt 2 m dx 



(3-1) 



When two particles initially at points x and x' are started with velocities v and v' and 
are driven by an identical sequence of random velocities v\ an at the same time intervals 
Tj, coordinates and velocities of them may accidentally draw closer to one another. The 
convergence of the two trajectories to the single final trajectory will depend on the evolution 
with a time of the small variances of the distance Ax, = x\ — Xi and velocity Avi = v[ — V{. 
Moreover, we investigate a transition from chaotic to nonchaotic behavior. Generally, such 
a transition may be detected from analysis of behavior of the neighboring trajectories and 
it is described by the Lyapunov characteristic exponents and KS metric entropy of the flow 
of trajectories in a given region of phase space [12, 13, 17-20]. 

From formal solutions x = x(xi,Vi,t) and v = v(xi,Vi,t) of equations (3.1) with initial 
conditions x = Xi and v — Vi at t — it follows an equation for Ax(t) and Av(t) at a time 
moment t: 

(twH^^Kt:) < 3 - 2 > 

where the matrix T is of the form 



I dx(xi,Vj,t) dx(xi,Vj,t) ^ 

dxi dvi 
dv(xi,Vi,t) dv(xi,Vi,t) 



dxi 



(3.3) 



Note, that the similar method of investigation is used in the theory of transition to chaos in 
classical systems [17, 18]. However, the motion in the form (3.2) and (3.3) is represented as 
the non-area-preserving tangent map, while classical dynamics of the conservative systems 
may be represented by the area-preserving maps. 

According to equations (3.1) and (3.3) matrix elements T xx and T xv satisfy the equation 



d 2 T x 
dt 2 



1 d 2 V 



m 



dx 2 



(3.4) 



while T vx = T xx , T vv = T xv and the initial conditions at t — are: 

T xx (^Xj / , v j, 0) T vv 1, T xv T vx 0, 

1 d 2 V 



T xx \Xi, Vi, 0) T vv 0, T xv 1, T v . 



m dx 2 



(3.5) 
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Therefore, the dynamics of the distance between the particles Ax and the difference of the 
velocity At> may be represented by the non-area-preserving mapping form of the equations 
of motion 

' Ax i+ i \ m/ s / Axi 



An 



i 1 )=T(a;x l ,^r l )(^). (3.6) 

In general, the intervals between the resets of the velocity 7* may be depending on the number 
of step i. 

Further analysis of the model may be based on the general theory of the dynamics of 
classical systems represented as maps [17-19]. Thus, for a = the Lyapunov exponent is 
defined as 

1 N 1 

\= \im —^2 — In T xx (xi,Vi,Ti) (3.7) 

00 i=l i 

and may be easily evaluated numerically. 

For a — 1 the map (3.6) is area-preserving and det T(l; Xi, Vi, Tj) = 1, while in general 
det T = a, TrT = T xx + aT vv and the eigenvalues /ii j2 of the T matrix are given by equation 

fi 2 - fiTrT + det T = 

which yields 



01,2 = 2 



T xx H - Q^T VV 

T \J (T xx + aT vv ) 2 — 4a 



So, the eigenvalues come in reciprocal pair, \i\\i2 = ol. For (T xx + aT vv ) 2 — 4a < the 
eigenvalues form a complex conjugate pair with = \/i 2 \ = yf®, otherwise the eigenvalues 
are real. 

Generally, the mapping T(a; t{) in (3.6) is depending on the starting coordinates X{ 
and Vi. Therefore, calculation of the mapping for n steps, T n = T(a; Xi+ n -i,Vi +n -i, r i+n _i) ■ 
T(a; Xj + „_2, fi+n-2, Ti+n-2) • • • T(a;Xi,Vi,Ti), and of the corresponding eigenvalues are 
complicated problems. Further we will evaluate the averaged quantities 

1 1^1 
<7i,2 = (-ln|/x 1>2 |) = lim — ^2 — In fj, lt2 {xi,Vi,Ti) (3.8) 

which are analogous of the averaged Lyapunov exponent (3.7), characterize the rate of the 
exponential increase of the separation of the two initially adjacent points and are related 
with the KS entropy of the system [17, 18]. Comparisons of the threshold values r c from 
the direct numerical simulations with those from the criterion 

^largest = (3.9) 

indicate to the usefulness of the quantities (3.8) for analysis of transition from nonchaotic 
to chaotic behavior and synchronization of the systems [21-23]. 

Therefore, theoretical analysis based on the mapping form of equations of motion for 
the distance between the particles and the difference of the velocity allows to simplify the 
problem of investigation of transition to nonchaotic behavior and results in the expressions 
for the criteria of the nonchaotic motion. Theoretical results agree well with the direct 
numerical simulations and indicate to the possibilities of generalization of the model, e.g. to 
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more degrees of freedom, for random values of the time intervals between the resets of the 
velocity and for systems driven by the random forces [21-23]. In paper [24] this method has 
been generalized and used for analysis of systems with repulsive force between particles, some 
scaling properties for the threshold reset time have been derived and it has been suggested 
that such convergence of chaotic orbits is a rather general phenomenon. 



4 Dynamics of quantum systems 
4.1 Two-level system 

Let's consider the simplest quantum dynamical process and the influence of frequent mea- 
surements on the outcome of the dynamics. Time evolution of the amplitudes a±(t) and 
a 2 (t) of the two-state wave function 

* = ai(t)*i + 02(*)* 2 (4.1) 

of the system in the resonance field (in the rotating wave approximation) or of the spin-half 
system in a constant magnetic field may be represented as 

f ai(t) — ai(0) cos |fii + ^(O) sin ^flt (A 2) 

[ a 2 (t) = iai(0)sin|fit + a 2 (0)cos|fit, ' ' ' 

where Q is the Rabi frequency. We introduce a matrix A representing time evolution during 
the time interval r (between time moments t = kr and t = (k + l)r with integer k) and 
rewrite Eq. (4.2) in the mapping form 



ai(k + 1)\ = A / a x (k) 
a 2 {k + l)J \a 2 {k) 



(4.3) 



where the evolution matrix A is given by 



( cosw isinw\ 1„ ,. 

A = . . r , if = -fir. (4.4) 



(4.5) 



Evidently, the evolution of the amplitudes from time t = to t = T = nr may be expressed 
as 

'ai(n) \ = An ( ai(0)' 
,a 2 (n)J \a 2 (0). 

One may calculate matrix A n by the method of diagonalization of the matrix A. The result 
naturally is 

A „ = / cos nip isinn(p\ ^ g . 

\i sin nip cos np ) ' 

Note that mp = \QT. 

Equations (4.2)-(4.6) represent time evolution of the system without the intermediate 
measurements in the time interval -j- T. If at t — the system was in the state i.e. 
ai(0) = 1 and a 2 (0) = 0, and if QT = n then at the time moment t = T we would certainly 
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find the system in the state ^ 2 , i-e. it would be |ai(T)| 2 = and |a 2 (T)| 2 = 1, a certain 
transition between the states. 

Let's consider now the dynamics of the system with the intermediate measurements every 
time interval r. Measurement of the system's state in the time moment t = kr projects the 
system into the state y&i with the probability pi(k) =\ a\(k) | 2 or into the state ty 2 with the 
probability P2{k) —\ a 2 (k) | 2 . After the measurement we know the probabilities pi(k) and 
p 2 (k) but we have no information about the phases a±(k) and a 2 (k) of the amplitudes 

ai (k) = \ ai {k)\ e iai(k \ a 2 {k) = \a 2 {k)\e ia2(k \ (4.7) 

i.e. the phases a±(k) and a 2 (k) after every act of the measurement are random [25]. This 
results in the equation for the probabilities 



p 2 (k + l)J ' \p 2 (k) 
where 

M=f C ° S 2 ^ Sin /) (4.9) 
V sin if cos if J 

is the evolution matrix for the probabilities. The full evolution from the initial time t — 
until t = T with the (n — 1) equidistant intermediate measurement is described by the 
equation 

,P2{n)J \P2{0)J 

The result of calculation of the matrix M n by the method of diagonalization of the matrix 
M is 

= l/l + cos"2v? l-cos"2v?\ , . 

2 V 1 - cos n 2ip 1 + cos n 2<p J ' 1 ' ; 

From Eqs. (4.10) and (4.11) we recover the quantum Zeno effect [25-27]: if initially the system 
is in the state than the result of the evolution until the time moment T = nr = n/il 
(after the n -pulse) with the (n — 1) intermediate measurement will be characterized by the 
probabilities pi(T) and p 2 (T) for finding the system in the states ^1 and ^2 respectively: 

1(1 + cos" 2^0 ~ |(1 + e~£) ~ 1 - g - 1, (4 12) 




1 — cos™ 2(f) ~ |(1 — e~^) ~ 'j- — > 0, n — > 00. 



We see that results of equations (4.10)-(4.12) represent the inhibition of the quantum dy- 
namics by measurements and coincide with those obtained by the density matrix technique 
[26, 27]. This also confirms correctness of the proposition that the act of the measurement 
may be represented as randomization of the amplitudes' phases. Further we will use this 
proposition and the same method for the analysis of the repeated measurement influence 
for the quantum dynamics of multilevel systems which classical counterparts exhibit chaos. 
We restrict ourselves to the strongly driven by a periodic force systems with one degree 
of freedom. The investigation is also based on the mapping equations of motion for such 
systems. 
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4.1.1 Dynamics of multilevel systems 

In general the classical equations of motion are nonintegrable and the Schrodinger equation 
for strongly driven systems may not be solved analytically. However, mapping forms of the 
classical and quantum equations of motion greatly facilitates the investigation of stochasticity 
and quantum-classical correspondence for the chaotic dynamics. From the standpoint of an 
understanding of the manifestation of the measurements for the dynamics of the multilevel 
systems the region of large quantum numbers is of greatest interest. Here we may use the 
quasiclassical approximation and convenient variables are the angle 9 and the action I. One 
of the simplest systems in which the dynamical chaos and its quantum localization may be 
observed is a system with one degree of freedom described by the unperturbed Hamiltonian 
Hq(I) and driven by periodic kicks. The full Hamiltonian H of the driven system may be 
represented as 

H(I,9,t) = H (I) + k cos 9 J2Ht~ 3 T ) ( 4 - 13 ) 

j 

where r and k are the period and strength of the perturbation, respectively. 

For the derivation of the quantum equations of motion we expand the state function 
i[>(9,t) of the system through the quasiclassical eigenfunctions, (p n (9) = e md / \f{2-n) , of the 
Hamiltonian H , 

1>{9, t) = (27T)- 1 / 2 J2 a n (t)i- n e- ine . (4.14) 

n 

Here the phase factor i~ n is introduced for the maximal simplification of the quantum map. 
Integrating the Schrodinger equation over the period r, we obtain the following maps for the 
amplitudes before the appropriate kicks 

a m (t j+1 ) =e- iHo ^ T Y t an{t j )J m -n{k), tj = jr (4.15) 

n 

where J m (k) is the Bessel function. 

The form (4.15) of the map for the quantum dynamics is rather common: similar maps 
may be derived for the monochromatic perturbations as well, e.g. for an atom in a mi- 
crowave field [28]. A particular case of map (4.15), corresponding to the model of a quantum 
rotator H = I 2 /2, has been comprehensively investigated with the aim of determining the 
relationship between classical and quantum chaos. It has been shown that under the onset of 
dynamical chaos at K = rk > K c = 0.9816, motion with respect to / is not bounded and it 
is of a diffusion nature in the classical case, while in the quantum description diffusion with 
respect to m is bounded, i.e. the diffusion ceases after some time and the state of the system 
localizes exponentially. This phenomenon turns out to be typical for models (4.15) with non- 
linear Hamiltonians Hq(I) and for other quantum systems. The quantum interference effect 
is essential for such dynamics and it results in the quantum evolution being quantitatively 
different from the classical motion. Quantum equations of motion, i.e. the Schrodinger equa- 
tion and the maps for the amplitudes, are linear equations with respect to the wave function 
and probability amplitudes. Therefore, the quantum interference effect manifests itself even 
for quantum dynamics of the systems, the classical counterparts of which are described by 
nonlinear equations; chaotic dynamics of the later exhibit a dynamical chaos. On the other 
hand, quantum equations of motion are very complex as well. Thus, the Schrodinger equa- 
tion is a partial differential equation with the coordinate and time dependent coefficients, 
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while the system of equations for the amplitudes is the infinite system of equations. More- 
over, for the nonlinear Hamiltonian H (m) the phases' increments, H (m)r, during the free 
motion between two kicks while reduced to the basic interval [0, 27r] are the pseudorandom 
quantities as a function of the state's quantum number m. This causes a very complicated 
and irregular quantum dynamics of the classically chaotic systems. We observe not only very 
large and apparently irregular fluctuations of probabilities of the states' occupation but also 
almost irregular fluctuations in time of the momentum dispersion [25, 28]. 

However, the quantum dynamics of such driven by the external periodic force systems 
is coherent and the evolution of the amplitudes a m (tj + i) in time is linear: they are defined 
by the linear map (4.15) with the time independent coefficients. The nonlinearity of the 
Hamiltonian H (I), being the reason of the classical chaos, causes the pseudorandom na- 
ture of the increments of the phases, H (ui)t, as a function of the state's number m (but 
constant in time). These increments of the phases remain the same for each kick. So, the 
dynamics of the amplitudes a m (t,- +1 ) = \a m (tj + i)\ e iam ^ tj+1 ^ and of their phases, a m (tj +1 ), 
is strongly deterministic and non-chaotic but very complicated and apparently irregular. 
For instance, the phases a m (tj + i) are phases of the complex amplitudes, a m (tj+i), which 
are linear combinations (4.15) of the complex amplitudes, a n (tj), before the preceding kick 
with the pseudorandom coefficients, e -' lH o( m ) r J m _ n {k). Nevertheless, the iterative equation 
(4.15) is a linear transformation with coefficients regular in time. That is why, we observe 
for such dynamics the quasiperiodic reversible in the time evolution [28] with the quantum 
localization of the pseudochaotic motion. 

In paper [28] it has been demonstrated that this peculiarity of the pseudochaotic quantum 
dynamics is indeed due to the pseudorandom nature of the phases, H (m)r, in Eq. (4.15) as 
a function of the eigenstate's quantum number m (but not of the kick's number j). Replacing 
the multipliers exp [— iH (m)r} in Eq. (4.15) by the expressions exp [— i2irg m ], where g m is a 
sequence of random numbers that are uniformly distributed in the interval [0, 1], we observe 
the quantum localization as well [28]. The essential point here is the independence of the 
phases H (m)r or 2ng m on the step of iteration j or time t. 

4.1.2 Influence of repetitive measurement on the quantum dynamics 

Each measurement of the system's state projects the state into one of the energy state <p m 
with the definite m. Therefore, if we make a measurement of the system after the kick j but 
before the next (j + 1) kick we will find the system in the states ip m with the appropriate 
probabilities p m (j) = \a m (tj)\ 2 . 

In the calculations of the system's dynamics the influence of the measurements may be 
taken into account through randomization of phases of the amplitudes after the measurement 
of the appropriate state's populations. The phases of amplitudes after the measurement arc 
completely random and uncorrelated with the phases before the measurement, after another 
measurements and with the phases of other measured or unmeasured states. Therefore, after 
the full measurement of the system after the kick j, all phases of the amplitudes a m (tj) are 
random. So, this full measurement of the system's state influences on the further dynamics 
of the system through the randomization of the phases of amplitudes. This fact may be 
expressed by replacement in Eqs. (4.15) of the amplitudes a m (tj + i) by the amplitudes 
e l ^ m ^' +1 ) a m (tj+i) with the random phases /3 m (tj+i). The essential point here is that the 
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phases f3 m (tj+i) are different, uncorrelated for the different measurements, i.e. for different 
time moments of the measurement tj + ±. This is the principal difference of the random 
phases f3 m (tj+i) from the phases iJo( m ) r in Eqs. (4.15) which are pseudorandom variables 
as functions of the eigenstate's quantum number m (but not of the time moment tj+i). 

Instead of representing the detailed quantum dynamics expressed as the evolution of all 
amplitudes in the expansion of the wave function (4.14) we may represent only dynamics 
of the momentum dispersion ((rrij — m ) 2 ) = E ( m ~ m o) 2 \o-m (tj)\ 2 where m is the initial 

m 

momentum (quantum number). Such a representation of the dynamics is simpler, more 
picturesque and more comfortable for comparison with the classical dynamics. 

Theoretically differences of dynamics without measurement and with the measurement 
may be understood from the iterative equations for the momentum dispersion: 



(m j+ i - m ) 2 ) = ( m - m o) \a m (t j+ i)\ , (4.16) 

m 

where 

\a m (tj+i)\ 2 = J ™-n (k) J m _ n ' (k) a n (tj) a*, (tj) . (4.17) 



Substitution of Eq. (4.17) into Eq. (4.16) yields 



(K+i - m ) 2 ) = E (m - m ) 2 J^„ n (k) \a n (tj) 12 



m,n 

\2 



+2 E E (m- mo) J m - n (k) J m _ w (k) Re a n (tj) a*, (tj) 

m.n i L 



n <n 



(4.18) 



For the random phase differences of the amplitudes a n (tj) and a\ (tj) with n ^ n, which 
is a case after the measurement of the system's state, the second term of Eq. (4.18) on the 
average equals zero. Then from Eq. (4.18) we have 



((m j+1 -m ) 2 ) = EK (tj)\ 2 J2(m-m ) 2 J 2 t _ n (k) 

n m 

= E \a m (tj)\ 2 (m 2 - m 2 + f ) = ((rrij - m ) 2 ) + f . 
In the derivation of Eq. (4.19) we have used the summation expressions 



(4.19) 



J2mJl- n (k) = and £ m 2 J 2 m _ n (k) = n 2 + ^ . (4.20) 

m m " 

Therefore, according to Eq. (4.19) for the uncorrelated phases of the amplitudes a n (tj) 
and a*, (tj) with n ^ n the dispersion of the momentum as a result of every kick increases 
on the average in the magnitude k 2 /2, the same as for the classical dynamics. Thus, we 
reveal that repetitive measurement of the multilevel systems with quantum suppression of 
classical chaos results in the derealization of the states superposition and restoration of the 
chaotic dynamics. Since this effect is reverse to the quantum Zeno effect we have called this 
phenomenon the 'quantum anti-Zeno effect 7 [25]. 
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5 Conclusions 



In this paper on the concrete examples we have demonstrated usefulness of the mapping 
equations for analysis of dynamics the nonlinear systems: transition to chaotic behavior in 
Hamiltonian systems, synchronization of chaotic systems driven by identical noise and effect 
of repetitive measurements for quantum dynamics. 
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